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I. Introduction 


The purpose of this report is to present the progress that has been made in the 
development of a model of manganese and iron fluxes from sediments. It is based on the 
sediment model (Di Toro and Fitzpatrick. 1993) developed for the Chesapeake Bay water 
quality model (Cerco and Cole, 1993). This report assumes that the reader has a 
familiarity with the concepts and formulations developed ther^n. 

This report concentrates on the development the manganese flux model. As shown 
below, the chemistry of iron and manganese are suffidentiy similar so that it is expected 
foat this model formulation can be applied to either. 

Two manganese data sets have been located, one of which is quite substantial. 

The equivalent data for iron, which indudes aerobic fluxes as well as anaerobic fluxes, 
does not appear to be available. Since the manganese data is much larger and more 
interesting it appears reasonable to use these sets for the initial development 



II. Chtmistry 


A. MangaiwM 

1. Spcctos and Radox Tianafbrmatfona 

The chemistry of manganese in natural waters and sediments has been studied for quite 
some time [Stumm and Morgan. 1970]. Manganese exists in two valances states: Mn(ll) 
in anoxic waters and Mn(iV) in oxic waters. Mn(IV) is very insoluble and forms manganese 
oxide, MnOjCs), which is the predominant form of manganese in oxic surtece waters. It 
usually exists as a coating on particles [Jenne.1968]. As the particles settle to the 
sediment, manganese is also transported providing a source of manganese to the 
sediments. In the oxic layer of the sediment. Mn02(s) is stable. However, particle mixing 
causes particles to be transported to the anaerobic layer of the sediment where 
manganese oxide is thermodynamicaliy unstable and a reduction reaction occurs. Mn(IV) 
is reduced to Mn(li). For this to occur, two electrons are required as can be seen from the 
following half reaction: 

MnOV)*** + 2e--►MnOO^^ (1) 

Thus the reaction requires an electron donor. The primary source of electrons in 
sediments is organic matter, CHjO, so that the reduction reaction can be written: 

Mn(IV)"* + ICHjO + IHjO -►Mn(ll)2* + IcO^ + 2H* (2) 

This reaction follows from the half reaction for oxygen: 

O 2 +4H*+4e*->HjO (3) 

and the oxygen equivalents of organic matter 
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CHjO + O 2 -^COj + HP 


(4) 


Since it is actuaiiy Mn02(s) tiiat is being reduced, the substitution of the equation 
corresponding to the formation of MnOjCs) yields: 

Mn(IV)"* + 2HjO -►MnOj(s) + 4H* (5) 


so that the complete reduction reaction of manganese oxide to Mn(ll) is: 

MnOj(s) + IcHjO + 2H* -^Mn^* + IcO, + 1 h,0 
2 2 ^ 2 ^ 


(6) 


In contrast to Mn(IV), Mn(ll) is more soluble and exists in the mgA. range in 
sediment pore waters. As a consequence it can diffuse to the oxic layer of the sediment 
where it is subject to oxidation. The oxidation of Mn(ll) to Mn(IV) occurs via the loss of two 
electrons: 

Mn(ll)2*^Mn(lvr+2e- (7) 

For oxygen as the electron acceptor, the overall reaction can be found using the half 
reaction for oxygen, eq.(3),: 

Mn(ll)2* +102 + 2H’ ->Mn(IV)** + HjO (8) 

followed by the precipitation of manganese oxide, eq.(5): 

Mn(ll)2* + + H 2 O -♦Mn02(8) + 2H* (9) 

This is the reaction that occurs in the aerobic layer. The kinetics of this reaction have bem 
examined [Morgan, 1967] and found to be slow in the normal pH ranges of surface waters. 
However the reaction can be bacterlally mediated and proceed more rapidly [Jaquet et al., 
1982]. 
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2. SotuWiity 


In oxic waters, Mn(IV) is very insoluble and manganese oxide, Mn02(s), is the 
predominant species. In anoxic waters, two solid phase species may exist: manganese 
carbonate, MnCOaCs) (rhodocrosite), and manganese sulfide, MnS(s). Since iron sulfide is 
also present in sediments, and it is more insoluble than MnS(s), it is unlikely that 
manganese sulfide is present We will assume below that the solubility of Mn(ii) is 
controlled by MnCOjCs). Thus not all of the Mn(ll) that is formed by the reduction of MnOj 
is in dissolved form. Some of it precipitates to form manganese carbonate. Therefore the 
transfer of Mn(ll) from the anoxic to the oxic layer occurs via particle mixing which 
transports MnC 03 (s) as weil as diffusion of soluble Mn(ll). 

B. Iron 

1. Soedes and Redox Transformations 

Like manganese, the chemistry of iron in natural waters and sediments has also 
been studied for quite some time [Stumm and Morgan, 1970]. Iron exists in two valances 
states: Fe(ll) in anoxic waters and Fe(lll) in oxic waters. Fe(lli) is very insoluble and forms 
iron oxyhydroxide, FeC>OH(s). Like manganese, it usually exists as a coating on particles. 
However, unlike manganese, there are other forms of iron that exist as particles in natural 
waters. Since the crust is approximately two percent iron, particles that runoff into r^atural 
waters contain a large amount of iron. As the particles settle to the sediment, iron is 
transported as well. This Is the source of iron to the sediments. 

Not all iron is reactive in sediments. It is convenient to denote the reactive portion 
of oxic iron as FeOOH(s) (Goethite), and to assume that it includes iron hydroxide as well 
since: 

Fe(OH) 3 (s) -►FeOOH(s) + HjO (10) 

The term iron oxyhydroxide is meant to denote the sum of FeOOH(s) and Fe(OH) 3 (s). 
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In the oxic layer of the sediment, FeOOH(s) is stable. However, as particle mixing 
causes particles to be transported into the anaerobic layer of the sediment, iron 
oxyhydroxide is thermodynamically unstable and a reduction reaction occurs. Fe(lll) is 
reduced to Fe(il). For this to occur, one electron is required as can be seen from the 
following half reaction: 

Fe(lll)®* + e-->Fe(ll)2* 

Thus it requires an electron donor. Again the source of electrons in sediments is 
organic matter, CHjO. Thus the reduction reaction can be written: 

Fe(lll)3* + IcHjO + IkO ->Fe(ll)2^ + IcO, + H* (12) 

4 4 4 

This reaction follows from eq.(3-4) and eq.(11). Since it is actually FeOOH(s) that is being 
reduced, the substitution of the equation corresponding to the formation of FeOOH(s): 

Fo(lll)3^ + 2HjO ->FeOOH(s) + 3H* (13) 

yields: 

FeOOH(s) + IcHjO + 2H^ -►FeOI)** + IcO^ + IHjO (14) 

4 4 4 

as the final reduction reaction of iron oxyhydroxide to Fe(ll). 

By contrast to Fe(iil), Fe(il) is more soluble and exists in the low mg/L range in 
sediment pore waters. As a consequence it can diffuse to the oxic layer of the sediment 
where it is subject to oxidation. The oxidation of Fe(ll) to Fe(iii) occurs via the loss of one 
electron which is the reverse of eq.(11). With oxygen as the electron acceptor, the overall 
reaction can be found using the half reaction for oxygen, eq.(3), 

Fe(ll)2* + lo. + H* -^Fe(lll)3’ + IHjO (15) 

4 2 

followed by the precipitation of iron oxyhydroxide, eq.(13) to yield the overall reaction: 
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( 16 ) 


Fe(ll)2* +IO 2 +2HjO ->FeCX)H(s) ♦2H* 

4 2 

This is the reaction that occurs in the aerobic iayer. The kinetics of this reaction have been 
examined and found to be slow in the normal pH ranges of surface waters. However, like 
manganese oxidation, this reaction can also be bacterially mediated and proceed more 
rapidly. 

Thus the chemistry of manganese and iron are quite similar. The oxidized forms 
are bofo insoluble and form oxides. The reduced forms are soluble in the mg/L range. 

Their concentrations in pore water are regulated by ~olid phases. Their flux to the 
overlying water is controlled by the extent that the reduced forms are oxidized in the 
aerobic layer, or escape as fluxes. Responses to lowered dissolved oxygen appear to be 
similar (Sunby et al., 1986). Hence it is expected that the formulations for manganese 
developed below can be applied to iron. 
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III. Partitioning Modal of Manganese Fluxes 


A. Model Formulation 

The model structure is shown in Fig.3.1. It is formuiated in terms of manganese but 
to equaliy appiies to iron since the mechanisms are anaiogous. There are four dependent 
variables: Mn(ll) and MnOjCs) in layers 1 and 2. These correspond to the total Mn(ll) and 
Mn(IV) in each layer. The source of manganese to the sediment is the settling of 
particulate manganese oxide with flux, Jwnoz. from \he overlying water. Two reactions occur 
in the aerobic layer. Mn(ll) partitions to form manganese carbonate. The reaction is 
parameterized with a linear partition coefRcient, x,. in addition, Mn(ii) is oxidized to 
Mn02(s) following eq.(9) with first order rate, 

Two reactions also occur in the anaerobic layer. Mn(ll) partitions to form 
manganese carbonate which is parameterized with a linear partition coefficient, Xj. This 
may be different from x^ in the aerobic layer due to the differences in pH and alkalinity. In 
addition, MnOjfs) is reduced to Mn(ll) following eq.(6) with first order rate, kwn. 2 - 

The mass transport between the overlying water and layer one is via the surface 
mass transfer coeffident, which is set equal to s = SOD/OjfO), the ratio of the 
sediment oxygen demand and the overlying water dissolved oxygen concentration, as in 
the previous modeis (Di Toro and Fitzpatrick, 1993). Partide mixing with mixing veiodty 
Wi 2 ^nd diffusive exchange with mass transfer coeffident {<^12 are induded as before, as is 
burial with sedimentation veiodty W 2 . 

B. Equations and Solutions 

The mass balance equations for the model follow from the reactions and transport 
processes discussed above. They are: 


7 





SEDIMENT 


MANGANESE FLUX MODEL 


WATER COLUMN 


Mn(ll)(0) 


SURFACE MASS TRANSFER: Klqi 
MANGANESE SOURCE: ^MnO^Cs) 

Itl 

PARTmONINQ: Mn(ll) ► MnCQsCs) 

*^Mn,1 

-► Mn02(6) 


OXIDATION: 


Mn(ll) 


PARTICLE MIXING 
Wl2 


DIFFUSION 

Kl12 


Cl 


PARTITIONING: Mn(ll) ^ ► MnC 03 (s) 

*^Mn,2 

REDUCTION: Mn02(8) -► Mn(IO 

SEDIMENTATION 
W2 


t 

BURIAL 
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Figure 3.1 









Layer 1 Mn(ll): 


0 - H, ^--^.11 - - 8 (f,,Mn(1) - Mn(0)) + (f«Mn(2) - f,,Mn(1)) 

dt 


+ w„(fpjMn(2) -fp,Mn(1)) -k„„.,f«Mn(1) 


(17) 


Layer 2 Mn(ll): 

0 --K,,2(f«Mn(2) -f,,Mn(1)) - (fp2Mn(2) -fp,Mn(1)) 

at 


-WjMn(2) +k^,„jMn02(2) 


(18) 


Layer 1 MnOji 

0 . H, - !i"^ .k.,,t.,Mn(1) *w„(MnO,(2)-MnO,(1)) 
cK 


(19) 


Layer 2 MnOj*. 


, dMn02(2) 
'^““dT” 


” *Silln .2 Mn02(2) 


-w„(Mn02(2) -Mn02(1)) -W2Mn02(2) 


( 20 ) 


where the particulate and dissolved fractions are computed from the partition coefficients 
and the concentration of solids, m, in eac^ layer 




1 

1 +mit 



( 21 ) 
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The solutions to these steady state equations are: 


Mn(2) » 


''MnO, 'Hla2 


^d1 ^liln.2 8^12 +f„s) +k„„_j) 


(22) 


where: 


ri2 “ 


^d2KLij +fp2W,2 


^di(KLi2 *^8) +fp,w« 


(23) 


The layer 1 solutions is: 


Mn(1) =r,2Mn(2) 
and the manganese flux is: 


(24) 


J[Mn] »s(f^,Mn(1)-Mn(0)) 


(25) 


This model will be compared to observations in the next section. 

C. Manganese Flux Data 

Two data sets wiii be examined using this model. This first is a reiativeiy smaii 
number of manganese and nutrient flux measurements made at three Long Island Sound 
stations (Aiier. 1980). The second is a iarge number of manganese (Hunt and Keiiy, 1988) 
and nutrient flux measurements made at the MERL mesocosms (Nixon et ai., 1986). The 
analysis technique is to examine the relationship between observed and computed 
ammonia, J[NHJ. and the manganese, J[Mn]. fluxes. As shown in Fig.3.2 there is a 
proportional relationship between these two fluxes. 
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J(Mii2I (mg Mn/m2-d) 


Manganese and Ammonia Flux 


Long Island Sound 



1000.0 


MERL Mesocosms 













1 . Rthitenship iMtwfMn J|NHJ and 800 

The parameter in the manganese model equations that depends on the ammonia 
flux is the surtace mass transfer coefficient s = SOD/O^fO). This parameter controls both 
the rate of mass transfer from the sediment to the overlying water and the depth of the 
aerobic layer, As s increases, one would expect that J[Mn] would increase since both 
the rate of surface mass transfer increases, and the depth of the aerobic zone decreases. 
The latter effect decreases the residence time in the aerobic layer arKt makes the oxidation 
of Mn(il) to Mn 02 (s) less rapid, enabling more Mn(ii} to escape to the overlying water. 

The dependency of s on ammonia flux occurs because the ammonia flux and SOD 
are related. This can be seen from the following equations. The steady state relationship 
between ammonia diagenesis, and ammonia flux, J[NH 4 ] is (Dl Toro and Fitzpatrick, 
1993): 


J[NHJ -4-^ (26) 

8* +K»IH4 

where the effect of the overlying water ammonia concentration is assumed to be small. At 
steady state, nitrogen diagenesis can be estimated using SOD and the Redfield ratio: 

SOD.ao2„4 (27) 


8 o 2 .n 58 the Redfield ratio between O 2 and nitrogen. Since: 
s * SOD/ 02 ( 0 ) . (28) 

s can be estimated from ammonia flux by solving these equations for s as a function of 
JINHJ. 
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This provides the necessary connection. Using eq8.(26-28), the relationship between 
ammonia diagenesis, J^. ammonia flux, J[NHJ and SOD is: 


'^N 


3 d, 


+ J[NHJ 



where; 

0, -I2<4,.J[NH/ * 27 J[NHJ 0,(0)“ 


(29) 


•a«JINHJic*^,O,(0)(4ai,,, JINHJ ♦27 i4,,,O,(0)' ) 

This approximation assumes that all cartx)n diagenesis eventually becomes SOD, 
i.e. that losses due to sulfide fluxes and burial are negligible. From an analysis of 
Chesapeake Bay fluxes, these losses amount to no more than 25%. The second source of 
error is due to the time lags that occur between the production of oxygen equivalents by 
carbon diagenesis and their eventuai oxidation. These are caused by the formation and 
the subsequent oxidation of FeS(8). From an analysis of steady state version of the 
sediment model K is known that the time tag effect causes an error of approximately a 
factor of two bebtraen the carbon diagenesis estimated from SOD assuming steady state, 
and the actual flux (Rg.5.4, Di Toro and Fitzpatrick 1993). Thus although these errors are 
not negligible, the approximations can be used so long as the magnitude of the error 
involved is understood. 

D. Comparison to Data 

The model computations are compared to the Long Island Sound data in Fig.3.3. 
The model parameters are listed in Table 1 and the computer program in the MATLAB 
langiM^ that produced the results is listed in Appendix I. The model predicts an 
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Manganese Linear Partitioning Model 




JINH4J (mf N/m^-d) JINH4I (mg N/m^-d) 


Manganese vs Ammonia Flux 



JiNH4l (mg N/m^-d) 
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Figura 3.3 





increasing flux of manganese as ammonia flux increases, Fig.(3.3a). As ammonia flux, and 
therefore s increases, less manganese is buried and. therefore, more escapes as a flux to 
the overlying water. The partition coefficient can be adjusted so that the observed 
particuiate and dissolved manganese concentrations in the anaerobic layer can be 
reproduced. The data are plotted arbitrarily at an ammonia flux of 10 mg N/m^-day. The 
model reproduces the concentrations at tt)e Ngher ammonia flux., Fig.(3.3c-d). However, it 
is dear that the model cannot reproduce the magnitude of the observed J[Mn] - JINHJ 
relationship, Fig(3.3e). The discrepancy - almost two orders of magnitude - is too large to 
be attributed to the error assodated with using s estimated from JENHJ. Therefore an 
important mechanism is missing in the linear partitioning model. 

E. Conclusions 

The linear partitioning model cannot reproduce the observed relationship between 
manganese and ammonia fluxes. At steady state, there are only two possible pathways for 
manganese: either it escapes as a flux, or it is buried. Therefore, in order for the model to 
reproduce the observations it is necessary that it predids a higher degree of burial at low 
ammonia fluxes and a larger manganese flux to the overlying water at high ammonia 
fluxes. 

Perhaps the problem is with the linear partitioning assumption. The fraction of 
manganese that is particulate in the anaerobic layer is determined by the solubility of 
MnC 03 (s). This is not a linear partitioning process. Rather it is controlled by a chemical 
equilibrium between manganese and carbonate ion concentrations. In order to model this 
process, it is necessary to model the processes that control the carbonate concentration in 
sediments. A model for this process is discussed in the next section. 
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Table 1 


Parameter Values for Linear Partitioning Model 


m. 

1.0 

kg/L 

m, 

1.0 

kg/L 

^12 

0.01 

m/d 


0.012 

m/d 

Wj 

0.25 

cm/yr 

*1 

100 

L/kg 


10.000 

Ukg 

^Mn.1 

1.0 

m/d 

I^Mn.2 

1.0 

m/d 

^^PP®NH41 

0.15 

m/d 

0,(0) 

5.0 

mg/L 

®02.N 

2.54*5.68 

mg Oj/mg N 

'^Mn02 

100 

mg Mn/m^-d 










































IV. Caiciuin - Alkalinity Flux Modal 


Moat freshwater and marine sediments contain a large concentration of calcium 
carbonate. Typical values are 10 -100 mg CaCOJg or 1 to 10% of the dry weight, so that 
there is as much inorganic carbon as there is organic carbon in sediments. This quantity of 
calcium carbonate provides a large buffer system fOr both the pH and the carbonate 
concentration in sediment pore water. The objective of the model formulated in the section 
is to reproduce the observations of pore water and solid phase concentrations of calcium 
and alkalinity. 

A. ChemMiy and SimpiHleations 

A modei of a chemical system is specified by the components in the model, and the 
species formed by the components (Morel. 1983). This applies to systems that are open 
(Oi Toro, 1976), i.e. subject to mass transport, as well as the normal closed system 
considered in typical chemical calculations. The calcium carbonate system is well 
urkterstood and is specified by three components: calcium, alkalinity, and total inorganic 
carbon. However, using three components results in an equation set that is too complex to 
be solved conveniently. Therefore a simplified set of components and equations are 
required. 

The equation that determines the solubility of calcium carbonate is: 


[CtfUCC^l -K,c«03 

The concentration of Ca^* can be approximated with the concentration of total dissolved 
catdum, Ca. This ignores the complexes of calcium with bicarbonate and ottier ligands. 
The concentration of carbonate can be approximated using the definition of alkalinity: 
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[AlkJ .[HCO^ +2[C0^1 +[OH1 - [HI 


(32) 


which In the pH range of sediment pore waters (pH = 7 - 8 ) and for large enough alkalinity. 
[Aik] >0.1 meg/L, is well approximated by; 


[Aik] -[HCO;] 


(33) 


The reaction between HCO 3 is: 


H* + col' o HCO3 


(34) 


and the mass action law is: 
[HC03 

[HUCOil * ^ 
on 

[Aik] 

(HKCOn 


so that 


[CX^l . . 

Thus the soiubility mass action equation becomes: 
ICaJAIk] 

KJHl N..C.C03 


(35) 


(36) 


(37) 


(38) 


or. 

'"'***• apparent solubility constant. Note that the effect of decreasing pH, 

which increases H*. increases the apparent solubility constant, as it should since [C 03 ^*j is 
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ICaHAIkl . 

Ci003 k,[hi.Kc^ 


(39) 


decreasing as pH decreases. 

This simplification reduces the number of components to two: Ca and Aik, and 
requires only that the pH be specified. This is a significant reduction in the complexity of 
the equations that need to be solved. 

For a simple closed system, the mass balance equations are: 


[Ca] + [CaAlk] = Ca^ 

(40) 

[Aik] +2 [CaAlk] ^Alk^ 

(41) 


The two in the alkalinity equation arises from the two equivalents of alkalinity in each mole 
of CaCOaCs). Ca^ and Alk^ are the total concentrations of calcium and alkalinity in the 
system. Substituting these equations into the mass action eq.(39) yields: 

(CSt-C aAlkXAIk^ -2CaAlk) =Kc^ (42) 

which is a quadratic equation that is solved for CaAlk. This is then substituted in eqs.(40- 
41) to determine Ca and Aik. 

B. Sediment Model Equations and Solutions 

The model formulation is shown in Fig.4.1. The dependent variables are: Alk(1) and 
Alk(2), the total alkalinity in layer 1 and layer 2 respectively; Ca(1) and Ca(2), foe total 
calcium in foe same sequence; and CaAik(2). the calcium carbonate in layer 2. The mass 
balance and mass action equations are as follows. 
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SEDIMENT 


CALCIUM FLUX MODEL 
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Figure 4w1 








Layer 1 Alkalinity: 

0 - s (Alk(O) - Alk(1)) * (Alk(2) - 2 CaAlk(2) - Alk(1)) (43) 

Layer 2 Alkalinity: 

0 “ -Kli 2 (Alk{2) - 2 CaAlk(2)-Alk(1)) - w, (2 CaAlk(2)) + (44) 

Layer 1 Calcium: 

0 » s (Ca(0) - Ca(1)) Ca(2) - CaAlk(2) - Ca(1)) (45) 

Layer 2 Calcium: 

0 “ -Kti 2 (Ca(2) - 2 CaAlk(2) - Ca(1)) - Wj (CaAlk(2)) (46) 

CaCO, Solubility 

(Ca(2) - CaAlk(2)) (Alk(2) - 2 CaAlk(2)) = (47) 

where: 

Kcai. = K^ccoa^CHI (48) 


is the apparent solubility product. The terms: Ca(2) - CaAik(2) and Aik(2) - 2 CaAlk(2) are 
the dissolved caldum and alkalinity in layer 2. 

The solution to these equations can be found as follows. The variables: Alk(1), 
Alk(2), Ca(1), and Ca(2) are eliminated from the five simultaneous equations. 
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The resulting equation for CaAlk(2) is a quadratic equation of the form: 


a CaAlk(2)* + b CaAlk(2) + c « 0 (^9) 

where: 

a-2(8 (50) 

b --(8 ♦K,„)((2Ca(0) *Alk(0))K,„8 ■.J„(8 ■•K,„))Wj (51) 


= - -Alk(0)Ca(0))Ki„8 -J„Ca(0)(8 »K^„)) (52) 

The computer code that implements these equations is given in Appendix 3. 

It is remarkable that the solution to these five simultaneous equations is reduced to 
the solution of a quadratic equation, as it is for the simple closed system. This suggests 
that so long as the chemistry can be simplified to the point that the equations for a closed 
system are solvable, then the open system equations can also be solved. 

C. Application to Long Island Sound 

The data for this application comes from observations of fluxes from three stations 
in Long Island Sound (^ler, 1980a, 1980b). The parameters required for the calcium - 
alkalinity flux model are the usual transport parameters as well as those specific to these 
components. These are listed in Table 2. The computer program that produced the results 
Is listed in Appendix II. 

In addition to the source of alkalinity from the overlying water there is also a 
sediment source of alkalinity. This is computed from an ammonia flux using the equations 
given above to obtain s and J^,. Redfield stoichiometry is used to compute the carbon 
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diagenesis, J^. It is assumed that all carbon diagenesis reduces suitate to sulfide following 
the formula: 

2 CHjO + SC^’ ->2 HCO 3 + HjS + 2Hp (53) 

Therefore each mole of organic carbon reacted produces one equivalent of bicarbonate 
alkalinity. 

The only chemical parameter required for the model is the apparent solubility of 
CaAlk. Fig.4.2 presents the pore water data for dissolved calcium, Ca, and alkalinity, Aik, 
from three Long Island Sound stations. The product; [Ca][Alk] is also presented. It ranges 
from approximately 20 (mM)^ to over 40 (mM)^. This sets the range for Kc^k- 

With the parameters established, the model is evaluated by specifying an ammonia 
flux and computing the resulting concentrations. The additional parameter values are 
presented in Table 2. 


Table 2 

Parameter Values for the Calcium Carbonate Model 


Ca(0) 


326 


mg Ca/L 


Alk(O) 


91.8 


mg CaCOj/L 


Kk... 


CaC03 


10-30 


mM^ 


1. Results 

Fig.4.3 presents the results for = 10-30 (mM)^. Alkalinity increases as 
ammonia flux increases, Fig.4.3a. For = 30 (mM)^ the solubility of CaAlk is exceeded 
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CaAlk (mM^) C* (mg Ca/L) 


Long Island Sound Pore Water 



0 I-^^-1 

74 75 76 77 

Year 


Alkalinity 



Apparent Solubility Product 
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for J[NHJ > 10 mg N/m^-d and CaC 03 (s) starte to form, Fig 4.3e. At tfiat point the calcium 
flux changes ffwn zero to negative, Fig.4.3b. which is the flux into the sediment that is 
required to support the burial of CaCO,. Also the flux of alkalinity from the sediment to the 
overlying water increases less rapidly, Fig.4.3a, since a portion is being buried. For smaller 
solubility products, =10-20 (mM)^, CaC 03 (s) forms over the entire range of J[NHJ 
investigated, Fig.4.3e, and therefore caidum flux is to the sediment, Fig.4.3b. 

Pore water concentrations are compared in Fig.4.3c-d. The data is plotted at J[NHJ 
= 10 mg N/m^-d for convenience only. The annual average ammonia flux is somewhat 
larger. Alkalinity and caidum are reproduced for = 30 (mM)^. Caidum concentrations 
as essentially equal to the overlying water concentrations whereas the alkalinity is larger, 
due to the additional source that results from sulfate reduction. 

The concentration of CaC 03 (s) is computed to increase from less than 10 mg 
CaC 03 /g (1% of dry weight) to over 100 mg (10%). Long island Sound 

sediments contain between 25 and 200 mg/g. The model reproduces the observations of 
the lower values of the solubility constant. This is in contrast to the pore water results. 
Perhaps there is an additional source of caidum carbonate to foe sediments, e.g, CaC 03 (s) 
from bivalve shells, which accounts for foe additional CaC 03 (s). 

Nevertheless, foe model is reasonably successfoi in reprodudng foe general 
features of foe data. Pore water caidum concentrations are predicted to be dose to foe 
overlying water value, 326 mg/L, whereas alkalinity is predicted to be larger than foe 
overlying value of 100 mgA.. Caidum carbonate concentrations are predicted to be in the 
range of 10 to 100 mg/g which is foe range of foe observations. These results are 
obtained using transport parameters that are calibrated from Chesapeake Bay sediment 
flux data, suggesting that these parameters are representative for these Long Island Sound 
sediments. In particular, foe sedimentation velodty for these sediments is quite dose to 
0.5 cm^r used for these calculations. 
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V. MangarwM - Caickim - Alkalinity Flux Modal 


A model is formulated in whi^ the pore water and solid phase concentration of 
mangarvese are controlled by the solubility of manganese carbonate. To this end, the 
caidum carbonate model is included as part of the formulation. The rest of the model 
ptfdlels the linear partitioning model presented in Section III. 

A. Chamislry and Simplifications 

The simplification of the manganese chemistry parallels that for calcium. The 
concentration of Mn(ll)^* is approximated wth the concentration of total dissolved 
manganese. Mn. This ignores the complexes of manganese with bicarbonate and other 
ligarfos. The concentration of carbonate can be approximated using alkalinity as 
formulated in section III. The mass action law for the solubility of MnCOjCs) is: 

[MnllAlk] (54) 

Pore water data from Chesapeake Bay sediments, Fig.5.1, (Bricker et al., 1977) and Long 
Island Sound, Fig.5.2, (Aller, 1980) are shown. MnAlk apparent solubility products are 
similar for both sets, ranging from 0.1 to 1.0 mM^. 

1. Equations and Solutions for MnCO,(s) and CaCO,(s) 

The mass balance and mass action equations for a dosed system are: 


(55) 

(56) 
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Alkalinity mass balance: 

AIICt -2CaAlk -2MnAlk -Alk^ »0 
Caidum mass balance: 

Ca^ -Ca^ -CaAlk =0 
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Mmganese mass balance: 

Mn^ - Mn^ - MnAlk - 0 (57) 

(^OgCs) Solubility: 

Ca^AIK. -K.,c«. (») 

MnCOgCs) Solubility 

Mn-AIK. (59) 

where Aik^. Aik^ are the total and dissolved alkalinity, Mn^ and Mn^ are the total and 
dissolved manganese, and Ca^ and Ca^ are the total and dissolved calcium. 

Eliminating all the independent variables except MnAlk yields a quadratic equation 


of the form: 

a MnAlk^ + bMnAlk + c = 0 (®0) 

with coefficients: 

b B 2 (2 Kip_y„Ai, Mrij + 4 Mn^^ - 2 Ca^ K^MnAm ^.MnAik) 

c - -2 ( 2 Mn? - 2 Ca, Mn, + Aik, Mn, - l<..HnA.c) (®3) 

The equation for CaAlk is also quadratic: 

aCaAlk^+bCaAlk+c *0 (64) 
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b »(2MnAlk - 2Caj -Alkr) 


( 66 ) 


c » -2 Ca^ MnAlk - + Alk^ Ca^ (®7) 

it is remarkabie that adding manganese as a component and MnC 03 (s) as a soiid phase 
does not materiaiiy compiicate the solution. It requires only the solutions of quadratic 
equations. 

The solution procedure is to solve the quadratic equation for MnAlk which yields two 
solutions. For each of these solutions, two solutions for CaAlk are found. For each of the 
four possible solutions the dissolved concentrations are computed. Then the solutions are 
checked for the following conditions. Are ail concentrations positive? Are there no 
oversaturated solids? If both of these conditions are true, then the feasible solution has 
been found. If no feasible solution is found, then either CaAlk and/or MnAlk are zero and 
the solution is undersaturated with either or both of these solids. The equations for these 
cases are given in the Appendix. 

2. Results 

/^n example computation is presented below. The computer program that produced 
foe results is listed in Appendix III. The catdum concentration is kept constant, Ca = 300 
mg/L = 7.5 mM which is approximately foe concentration in overlying water for foe Long 
Island Sound sediment data. The solubility products are; (mM)^ and K^, MnAk ~ 

0.4 (mM)^. Total alkalinity is varied from 1 to 20 mM and total manganese from 0.05 - 0.5 
mM. These ranges are characteristic. Note that foe alkalinity is much larger than foe 
manganese concentrations. The tables below have alkalinity varying in foe y direction and 
manganese in foe x direction. 





The computed solid 

phase 

concentrations. MnAlk and CaAlk, are: 

MnAlk (mM) 


AlCr(mM) 

1.0000 

0 

0 

0 

0 

0 

0.699 

1.8206 

0 

0 

0 

0 

0.0309 

0.2094 

3.3145 

0 

0 

0 

0 

0.0777 

0.2566 

6.0342 

0 

0 

0 

0.0389 

0.1525 

0.3325 

10.9856 

0 

0.0045 

0.0504 

0.1230 

0.2382 

0.4207 

20.0000 

0.0198 

0.0490 

0.0952 

0.1685 

0.2847 

0.4688 

MrirCmM) 

0.500 

0.0792 

0.1256 

CaAlk (mM) 

0.1991 

0.3155 

0.5000 

AlCrCmM) 

1.0000 

0 

0 

0 

0 

0 


1.8206 

0.4101 

0.4101 

0.4101 

0.4101 

0.3844 

0.2348 

3.3145 

1.6152 

1.6152 

1.6152 

1.6152 

1.5548 

1.4146 

6.0342 

3.5210 

3.5210 

3.5210 

3.4971 

3.4266 

3.3132 

10.8856 

5.6321 

5.6309 

5.6190 

5.5998 

5.5689 

5.5184 

20.0000 

6.7445 

6.7429 

6.7404 

6.7364 

6.7299 

6.7195 

MivtmM) 

0.0500 

0.0792 

0.1256 

0.1991 

0.3155 

0.5000 


CaAlk forms as Alk^ is increased. Its concentration is not affected by the concentration of 
Mny. MnAlk also forms as Alkj and Mnr are increased. Its concentration is affected by 
both independent variables. 
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The solubility products are listed below. Note that the aqueous phase is 
undersaturated where no solid phase forms. 


AHir (mM) 

1.0000 

0.0500 

0.0792 

1.8206 

0.0705 

0.1118 

3.3145 

0.0850 

0.1347 

6.0342 

0.1257 

0.1992 

10.9856 

0.2677 

0.4000 

20.0000 

0.4000 

0.4000 

Mrir (mM) 

0.0500 

0.0792 


ANCr (mM) 

1.0000 

7.5000 

7.5000 

1.8206 

10.0000 

10.0000 

3.3145 

10.0000 

10.0000 

6.0342 

10.0000 

10.0000 

10.9856 

10.0000 

10.0000 

20.0000 

10.0000 

10.0000 

Miv (mM) 

0.500 

0.0792 


[Mn][A)k] (mM)' 

! 


0.1256 

0.1991 

0.3155 

0.1771 

0.2808 

0.4000 

0.2134 

0.3382 

0.4000 

0.3156 

0.4000 

0.4000 

0.4000 

0.4000 

0.4000 

0.4000 

0.4000 

0.4000 

0.1256 

0.1991 

0.3155 

[Ca][Alk] (mM)' 

2 


7.5000 

7.5000 

7.5000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

0.1256 

0.1991 

0.3155 


0.4000 

0.4000 

0.4000 

0.4000 

0.4000 

0.4000 

0.5000 


6.9756 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

0.5000 
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The dissolved concentrations are listed below. They are ail positive as required by 
the feasibility conditions. 


Dissolved Mn (mM) 


Alkr(mM) 


1.000 

0.0500 

0.0792 

0.1256 

0.1991 

0.3155 

0.4301 

1.8206 

0.0500 

0.0792 

0.1256 

0.1991 

0.2846 

0.2906 

3.3145 

0.0500 

0.0792 

0.1256 

0.1991 

0.2378 

0.2434 

6.0342 

0.0500 

0.0792 

0.1256 

0.1601 

0.1629 

0.1675 

10.9856 

0.0500 

0.0748 

0.0752 

0.0760 

0.772 

0.0793 

20.0000 

0.0302 

0.303 

0.304 

0.0305 

0.0308 

0.0312 

MnrtmM) 

0.0500 

0.0792 

0.1256 

Dissolved Aik 

0.1991 

(mM) 

0.3155 

0.5000 

AilCr(mM) 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.9301 

1.8206 

1.4105 

1.4105 

1.4105 

1.4105 

1.4105 

1.3764 

3.3145 

1.6993 

1.6993 

1.6993 

1.6993 

1.6820 

1.6433 

6.0342 

2.5132 

2.5132 

2.5132 

2.4982 

2.4550 

2.3885 

10.9856 

5.3535 

5.3502 

5.3163 

5.2627 

5.1784 

5.0465 

20.0000 

13.2358 

13.2081 

13.1644 

13.0951 

12.9854 

12.8118 

Mrif (mM) 

0.0500 

0.0792 

0.1256 

Dissolved Ca 

0.1991 

(mM) 

0.3155 

0.5000 

Akr(mM) 

1.0000 

7.5000 

7.5000 

7.5000 

7.5000 

7.5000 

7.5000 

1.8206 

7.0899 

7.0899 

7.0899 

7.0899 

7.1156 

7.2652 

3.3145 

5.8848 

5.8848 

5.8848 

5.8848 

5.9452 

6.0854 

6.0342 

3.9790 

3.9790 

3.9790 

4,0029 

4.0734 

4.1868 

10.9856 

1.8679 

1.8691 

1.8810 

1.9002 

1.9311 

1.9816 

20.0000 

0.7555 

0.7571 

0.7596 

0.7636 

0.7701 

0.7805 

MfirtmM) 

0.0500 

0.0792 

0.1256 

0.1991 

0.3155 

0.5000 
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B. S«difn«nt Model Equations and Solutiona 


Vhe structure of the manganese flux model is illustrated in Fig.5.3. The oxidation of 
Mn(ll) to MnOjCs) ccurs in the aerobic layer, and the reduction or Mn02(s) to Mn(ll) occurs 
in the anaerobic layer. The partitioning of Mn(li) in the anaerobic layer is controlled by the 
solubility of MnCOjCs). In order to calculate the carbonate concentration, the influence of 
CaCOjCs) must also be considered. Both these reactions are shown as the solubility of 
MnAlk and CaAlk, as discussed above. 

Hence, the equations for the sediment model are a combination of the equations for 
the linear partitioning model for manganese, with the addition of the equations for calcium, 
from the calcium - alkalinity model, and a solubility equation for MnAlk. They are listed 
below. 


1. Mass Balance Equations: 

Layer 1 Alkalinity; 

-s (Alk(1) - Alk(O)) + (Alk(2) - 2 CaAlk(2) - 2 MnAlk(2) - Alk(1)) = 0 (68) 

Layer 2 Alkalinity; 

-KLi2(Alk(2) -2CaAlk(2) -2MnAlk{2) -Alk(1)) 

-W 2 (2 CaAlk(2) + 2 MnAlk(2)) + = 0 


Layer 1 Calcium; 

-s Ca{1) + Kli 2 (Ca(2) - CaAlk{2) - Ca(1)) - s Ca(0) = 0 (70) 
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F 


Layer 2 Calqum: 

(Ca(2) - CaAlk(2) - Ca(1)) - CaAlk(2) =0 (71) 

Layer 1 Mn(ll): 

(-8(fj, Mn(1) - Mn(0)) + (Mn(2) - MnAlk(2) - Mn(1)) 

+ w„(f^MnAlk(2) -fpiMn(1)) Mn(1)) =0 (72) 

Layer 2 Mn(ll): 

( -Kl,j (Mn(2) - MnAlk(2) - Mn(1)) 

- (fa„ MnAlk(2) - fp, Mn(1)) - MnAlk(2) + k^^j Mn02(2)) = 0 (73) 

Layer 1 Mn02(s): 

Mn(1) + w,2 (Mn02(2) - MnOjCI)) + J„„o2 = 0 

Layer 2 Mn02(s); 

-k*,„ 2 Mn02(2) - w,2 (Mn02(2) - Mn02(1)) - Wj Mn02(2) = 0 (75) 

2. Mass Action Equations: 

Layer 2 CaC 03 (s) solubility: 

(Ca{2) - CaAlk(2)) (Alk(2) - 2 CaAlk(2) - 2 MnAlk(2)) - = 0 (76) 
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Layer 2 MnCOgCs) solubility: 

(Mn(2) - MnAlk(2)) (Alk(2) - 2 CaAlk(2) - 2 MnAlk(2)) - = 0 (77) 

3. Solution 

The first step in the solution is to solve for ail the dependent variables except 
CaAlk(2) and MnAlk(2) using the mass balance equations. The solutions are: 

Ai.,M - < Alk(2) - 2 CaAlk(2) - 2 MnAlk(2) + Alk(O) s) 

(K,„ -^s) 

Ca(1) = (-(CaAlk(2)K,,,) -^Ca(2)K,,, -KCa(O)s) 

(Kli 2 +s) 

Mn(1) » (-(Kli 2 MnAlk(2)) + Kl ^2 Mn{2) + Mn(0) s 

+ f,^ MnAlk(2) f,, K,,, + f,, kMnI + f,, s + f,, w„) (80) 
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Mn02(1) ■ ( ( kniij ^12 ^2 ) ^ *^Mn02 ^Mnl ■*■ *^Mn02 ® 

+ Mn(0) s MnAlk(2) Wj ) / (w„ (Is,^ s + k„„, Wj + 8 Wj))) (81 ) 

Mn02(2) = ( ( Jyn02 ^Mnl * '^Mn02 ® ■*" ^Mnl ^^^(0) S 

-ky„i MnAlk(2) W 2 ) / (ky^ s + ky„, W 2 + s W 2 )) (82) 

It is the layer two solutions that are needed for the mass action equations. These are 
written in the following form in order to substitute for various parameter groups and to 
isolate the dependency on CaAlk and MnAlk: 

Alk(2) *ja0 <«-jacCaAlk(2) jam MnAlk(2) (®^) 

Ca(2) = jcO + jccCaAlk(2) + Jem MnAlk(2) (®^) 

Mn(2) =JmO jmcCaAlk(2) jmm-*MnAlk(2) (®5) 

where the notation denotes the Jacobian (j) of the equation (a,c,m for Aik, Ca. Mn) with 
respect to nothing, (0, the constant term) or the variables (c, m for CaAlk, MnAlk). 





The Jacobians for alkalinity are: 


jaO =Alk(0) +jalk / KI12 +jalk / s 

jac «2 -2Wj / - 2 w 2 / 8 (®7) 

jam *2 -2Wj / Kli 2 -2^^! s (®®) 

For calcium they are: 

jc0*Ca(0) W 

icc^-w2/KI12-w2/s 

jcm^ (91) 
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and for manganese they are: 


jmO ■ ( ( ^MKXt ^12 Kyln2 til '^Mn02 ^n2 


■*■ til ^^MnOa ^n2 * ■*■ fjl ^12 ^n2 Mn( 0 ) 8 


^ ^dl ^Mnl ^Mb2 Mn(0) 8 + fp, Jmmo2 Klitf ^*^12 


^ikM„ 2 Mn( 0 ) 8 w,j +f,,K,„Mn( 0 ) 8 Wj + fp, Mn( 0 ) 8 w,j Wj) 


f (tdi f^i 2 (^Mni ^ Kinl ^*^2 * 


jmc *0 


jmm ® ( ( ^n2 ® ” tit* K«n2 ® ^12 '*’ ^d1 ^12 ^n1 ^2 ~ ^d1 ^12 ^n2 ^2 


“^d1 Kdnl Kdn2 '*'^2 ■*■ ^d1 ^12 * - fp, 8 Wj - fpj, fp, ky„, W,2 Wj 


■^dld ^d1 ® '*'^12 '*'^2 ” ^d1 ^12 ^ ~ td1 ^MbI ^ “ ^dl ® ~ ^pl ^^12 ) 
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f (^dA«(»W® +*Wj))) (94) 

An important raiationsNp among these coefficients is: 

jac ■2jcc (®®) 

jam »2Jcc (96) 

Using these equations and substituting into the mass action equations yieids 
quaMic equations for both CaAlk and MnAik. The form is: 

aCaAlk* ♦bCaAlk+c-0 (97) 

and similarly for MnAik. The coefficients for MnAlk(2) are denoted by the subscript Mn: 

"2k,po^ +2jmm^k,pQ^ '*‘i^*'^(*4k,pciAic “^k^^,^) 

+ 2k^^ +jcc(-2k^^ +2jmmk^^) (98) 

b^ -im0(-4kc^ +4jmmkc^ -2k„„** 

+ic0(2k„„«. -Zjmmk.,^) +ia0(-k„^ +jmmk„^ (99) 
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--(-2jmO*ICc^ *im0(-(i30k^^J +2jc0k„^)) 


( 100 ) 


and for CaAlk(2) which are denoted by the subscript Ca: 

a^ -2 -4jcc +2jcc? (101) 

be «ja0(-1 +jcc) +(-2 +2|cc)jc0 +2MnAlk(2) 


-4jccMnAfk(2) •^2jC(^MnAik(2) 


( 102 ) 


Ce - -( -(JaO jeO) + keA„ + jeO (2 MnAlk(2) - 2 jee MnAlk(2)) (103) 


The solution procedure is to solve the quadratic equation for MnAlk(2) which yields two 
solutions. For each of these solutions, two solutions for CaAlk(2) are found. For each of 
the four possible solutions the rest of the variables are computed. Then the solutions are 
checked for the following conditions. Are all concentrations positive? Are there no 
oversaturated solids? If both of these conditons are true, then the feasible solution has 
been found. If no feasible solution is found, then either CaAlk and/or MnAlk are zero and 
the solution is undersaturated with either or bofo of these solids. 

Once again, the solutions are found from the roots of quadratic equations. They are 
very similar to the batch reactor equations, although the coefficients of the equations are 
much more involved. Nevertheless, no essential complexity is introduced and the 
numerical solutions are straightforward. 

C. ResuNs 
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A pralimin«y calibration to the Long IsUhxI Sound data is sho^ in Fig.5.4. The 
adcKtional parameter values we presented in Table 3 and the computer program that 
produced the results is listed in Appendbc IV. 


Tables 

Parameter Values for foe MnCO, and CaCO, Model 



0.03 

m/d 

^Mn.2 

10 

m/d 

*1 

300 

L/kg 

^C«003 

20 

mM* 

^.MnOOS 

0.5 

mM* 


The observed pore water alkalinity, calcium, and CaCOjCs) are reproduced by foe model 
calculations. The data are plotted at J[NH 4 ]» 10 mg/m’-d for convenience. The fit is 
actually slightly better than foe model for CaCOjCs) alone, Fig.4.3. 

The manganese results are shown la Fig.5.5. Pore water manganese is reproduced 
reasonably closely (C) while foe computed concentration of MnCOjCs) is higher than foe 
observation (D). The manganese flux variation with respect to ammonia flux (A,B) is more 
prorxxjnced then foe linear partitioning model, Rg.3.3. However it is still not as 
pronounced as foe observations (E). 

D. Conclusions 
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Figur* 5.5 







The inclusion of carbonate solubility controls for manganese partitioning In the 
artaerobic layer has slightly improved the model results. However, there is still a 
substantial variation of J[Mn] with respect to J[NHJ that is unexplained. Perhaps it is a 
time variable effe^. This possibility is examined in the next section. 
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VI. Application to MERL Data 


A. Raviaw or MERL Data 

The MERL mesocosms are large outdoor tanks approximately 1.8 meters in 
diameter and have a water depth of 5.0 meters. Narragansett Bay water flows through at a 
rate to establish a detention time of approximately 30 days. The tanks have a mixer to 
control stratification. Sediments are obtained using a large box core which maintains the 
vertical orientation and the top 40 cm are placed into containers in the bottom of the tanks 
(Nbcon et al., 1986). 

The data to be analyzed LwIow comes from the Nutrient Addition Experiment. Its 
purpose was to examine the consequences of nutrient enrichment to coastal estuaries. The 
duration was approximately 2 1/4 years over three calendar years. The nutrient dosing 
was increased in a geometric series, IX, 2X, 4X, 8X, 16X, and 32X, in addition to three 
control tanks. The nutrients added were inorganic nitrogen, phosphorus, and silica in a 
molar ratio of 12.8 N : 1.0 P : 0.91 Si to match the stoichiometry of sewage entering the 
bay (Nixon et al., 1986). Areal loading rates of total nitrogen to the tanks varied from 23 
mg N/m’-d for the controls, 63 mg N/m^-d for IX, 103 mg N/m^-d for 2X and so on 
geometrically to 1308 mg N/m’-d for the 32X (Kelly et al., 1985). As a result, mean annual 
water column DIN increased from 56 to 4200 ug N/L; mean annual Chia ranged from 4 to 
70 ug/L, and total system carbon production ranged from 0.55 to 2.2 g C/m^-d (Nixon et al., 
1986). Note that the areal loading rate increased approximately 57 fold, the DiN 
concentration by approximately the same ratio (75), whereas the chlorophyll a inaeased by 
only 17 fold and the total system carbon production increased by only 4 fold. 

Sediment processes were also examined during the experiment. Sediment oxygen 
and nutrient fluxes, pore water and solid phase concentrations were also measured. In 
addition, manganese flux and sediment compositional data were also collected (Hunt and 
Kelly, 1988). These data are examined below. 
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B. Application ofttw Nutrient Flux Model 


The stand alone version of the Chesapeake Bay Sediment Model was applied to the 
MERL sediment data. In the stand alone version of the model, the average annual 
depositional flux of particulate organic nitrogen to the sediment, JpoN. is specified externally. 
Values are chosen to fit the measured ammonia fluxes. No kinetic or transport parameter 
values were changed initially. The required depositional fluxes are shown in Fig.6.1. 

Fluxes vary from less than 50 to 130 mg N/m^-d, less than a 3 fold variation. 

The model is run as follows. The sediments are initially equilibrated using the 
depositional flux that characterizes the controls during the experiment. That is, it is 
assumed that the state of the sediments in ttie year when they were collected were similar 
to that in the subsequent years as indicated by the fluxes in the control mesocosms. Thus, 
initially all the sediments have the same concentrations as the controls. The differences 
occur due to the increased loading during the experiment. 

The resulting ammonia fluxes are shown in Fig.6.2. The seasonal variation is 
reasonably well reproduced except for 32X. However, the negative fluxes in the winter, 
which increase as loading increases, are not captured. 

The oxygen fluxes are shown in Fig.6.3 The results are remarkably good. The absolute 
magnitudes and the seasonal variations are reproduced up to 32X. Since the manganese 
model requires the surface mass transfer coefficient: s = SOD/O^iO), it is important that the 
oxygen flux be well reproduced by the model. 

C. Application <rfthe Manganese Linear Partitioning Model 

The linear partitioning model is described in Chapter III. A time variable version of 
the model has been implemented which utilizes the same equations as the steady state 
model. The same implicit integration technique is used for manganese as for the nutrients. 
It requires the solution of four simultaneous linear equations instead of the two for the other 
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variables, since the equations for Mn(il) and MnOjCs) are coupled. This presents no 
difficulty, however, and the implementation is straightforward. 

The model is run as follows. The nutrient portion is run as described above. For 
manganese, the sediments are initially equilibrated using a constant depositional flux. 

Then this flux is continued throughout the three years of the experiment. Thus, it is 
assumed that water column processes during ttie experiment do not change the 
depositional flux of manganese to the sediment. This is oniy an approximation since the 
overlying water concentrations vary as dosing varies. Fig.6.4 presents averaged seasonal 
data. The more highly dosed tanks usually have larger concentrations than the controls 
with the exception of 32X. This is most apparent in the earlier period of the experiment. 
However, overall the variation is not too lar^e so that the assumption of a constant flux 
from the overlying water appears to be a reasonable first approximation. 

A comparison of the model and observed fluxes is shown in Fig.6.5. The model is 
incapable of producing the strong seasonal variations in fluxes that are observed, 
especially at the higher loading rates. This is consistent with the findings of tfie steady 
state model. This is clearly illustrated In Fig.6.6 which compares the observed relationship 
between observed (left hand side) and modeled (right hand side) manganese flux and 
various other variables. The relationship to ammonia flux (bottom) demonstrates that the 
model does not reproduce the magnitude of ttie variation that is observed. 

The pore water and solid phase concentrations are compared in Fig.6.7. The pore 
water concentrations are reproduced quite nicely. The sediment concentrations are 
approximated by the minimum annual concenb'ations. Since the model generates an 
increasing manganese flux as the dosing increases it calculates a larger loss of sediment 
manganese for the more highly dosed tanks. The reason this occurs is that the higher 
dosed tanks have a larger SOD and, therefore, a larger mass transfer coefficient and 
thinner aerobic layer. Thus less Mn(ll) is oxidized and more escapes from the sediment. 
Since the quantity that is being deposited is constant across the dosing gradient, the higher 
dosed tanks lose more manganese. However the calculated effect is not as large as is 
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Figure 6.4 
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FIgur* 6.6 
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ob8«rvect since the model cannot reproduce the magnitude of the variation that is observed 
between increased loading and manganese fluxes. 

The strong relationship between manganese and ammonia flux is striking because 
ammonia is generated by the decomposition of organic matter, whereas manganese is not 
How, then, does this strong covariation come about? The model already accounts for the 
variations that would be e^q^ected from the varying thickness of the aeroixc layer. There 
must be another mechanism that is somehow related to diagenesis. 

Consider the following possibility. As diagenesis occurs organic matter decompK>ses 
and the pH of pore water drops. The extent depends on the exactly which reactions are 
assumed to occur. A drop in pH corresponds to an increase in the solubility of MnAlk, see 
eq.(S4). This would cause less manganese to be buried and more to escape. It is not 
krK)wn whether this mechanism can explain the observed seasonal variations. However, it 
appears to be the next logical step in the development of the manganese flux model. 


D. Conclusions and Recommendations 

At the present stage of model development, the manganese model captures some 
of the features of the data sets examined. However, it is incapable of reproducing the 
magnitude of the seasonal variation of manganese fluxes. This is true whether it is 
examined as a function of ammonia flux, as in the Long Island Sound data, or as a time 
series, as in the MERL data. Since the model Incorporates changing thickness of the 
aerobic layer, which is the usual explanation of the seasonal variation, there must be 
another mechanism operating. The next step in foe model development is to implement 
foe mechanism that affect manganese carbonate solubility - decreasing pH with increasing 
diagenesis - in order to reproduce foe observed seasonal variation. 
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